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ABSTRACT 


A  model  forest  canopy  was  designed  to  simulate  the 
meteorological  characteristics  of  typical  live  forests. 
Measurements  were  made  of  velocity,  turbulence,  drag,  and 
gaseous  plume  behavior.  Flow  properties  are  compared  with 
recent  field  measurements.  Ground  penetration  in  the 
initial  fetch  region  results  in  strikingly  different  stream¬ 
line  motion  as  compared  to  wind  motions  within  the  equilib¬ 
rium  regions.  Measured  values  of  the  vertical  eddy  dif¬ 
fusion  coefficient  are  shown  to  predict  plume  behavior  in 
tie  equilibrium  region  very  well  if  a  correction  is  in- 

j £ 

eluded  for  the  ratio  >  1.0 

Kz 

Ventilation  of  an  elevated  line  source  into  the  canopy 
region  is  compared  with  a  simple  one-dimensional  model. 


TABLE  OF  CONTENTS 


Page 

ABSTRACT .  ii 

LIST  OF  FIGURES .  iv 

LIST  OF  TABLES .  vi 

1.  INTRODUCTION  .  1 

2.  MODELING  OF  A  FOREST  CANOPY  .  4 

3.  EXPERIMENTAL  EQUIPMENT  AND  PROCEDURES  .  7 

3.1  Wind  Tunnel  and  Canopy  Arrangement  ....  7 

3.2  Velocity  and  Turbulence  Measurements  .  -  .  7 

3.3  Concentration  Measurements  -  Helium 

Tracer  Gas .  8 

3.4  Concentration  Measurements  -  Kr-85 

Tracer  Gas .  9 

4.  EXPERIMENTAL  RESULTS .  11 

4.1  Typical  Velocity  and  Turbulent  Intensity 

Profile  Results  .....  .  11 

4.2  Diffusion  Plums  Results  .  13 

4.3  Eddy  Diffusion  Coefficient .  16 

4.4  Forest  Penetration  Models  .  19 

5.  CONCLUSIONS .  23 

REFERENCES .  25 

TABLES .  29 

FIGURES . . .  36 


LIST  OF  FIGURES 

Figure  Page 

1  Wind  tunnel  arrangement  .  36 

2  Model  plastic  forest  .  37 

.  3  Helium  detectio:  r  . 38 

4a  Krypton-85  detection  system  -  source  ...  39 

4b  Krypton- 8 5  detection  system  -  detector  .  .  40 

5  Drag  coefficient  of  live  and  model  trees  .  41 

6  Wake  characteristics  of  live  and  model 

trees . 42 

7  Shear  plate  drag  for  model  forest  canopy  ,  43 

8  Velocity  profiles  .  41 

9  Turbulence  intensities  .  45 

10  Comparisons  with  winter  forests  .  46 

11  Comparisions  with  summer  forests  .  47 

12  Velocity  defect  comparison  .  48 

13  Diffusion  -  Isoconcentration  profiles 
zs  =  0.0  cm 

xs  =  0.0  m . 49 

14  Diffusion  -  Isoconcentration  profiles 
zs  --  10.0  cm 

xs  =  0.0  m . 50 

15  Diffusion  -  Isoconcentration  profiles 
zs  =  0.0  cm 

xs  =  6.0  m . 51 

* 

16a  Diffusion  -  Isoccncentration  profiles 

zs  =  10.0  cm . 52 

16b  xs  =  6.0  m . 53 

17  Diffusion  -  Isoconcentration  profiles 
zs  =  18.0  cm 

xs  =  6.0  m . 54 

18  Diffusion  -  Isoconcentration  profiles 
zs  *  27  cm 
xs  =  6.0  m 


iv 


55 


LIST  OF  FIGURES  (continued) 


Figure 

Page 

19 

Ground  concentration  vs  downstream 
distance  . 

56 

20 

Dimensionless  above  canopy  concentration 
profiles  . 

57 

21 

Eddy  diffusion  coefficient  -  mass  .... 

58 

22 

Below  canopy  dimensionless  eddy  diffusion 
coefficient  profiles  . 

59 

23 

Analytical  check  on  ground  concentration 
variation  . 

60 

24 

Cross-section-isoconcentration  profiles  . 

61 

25 

One  dimensional  penetration  model  .... 

62 

26 

Coefficient  m  rs  Environmental  Index  .  . 

63 

v 


LIST  OF  TABLES 


Table  Page 

1  Concentration  profiles  of  diffusion 

in  the  plastic  tree  canopy . 29 

2  Concentration  profiles  of  diffusion 

in  the  plastic  tree  canopy . 30 

3  Concentration  profiles  of  diffusion 

in  the  plastic  tree  canopy . 31 

4  Concentration  profiles  of  diffusion 

in  the  plastic  tree  canopy . 32 

5  Concentration  profiles  of  diffusion 

in  the  plastic  tree  canopy . 33 

6  Concentration  profiles  of  diffusion 

in  the  plastic  tree  canopy . 34 

7  Concentration  profiles  of  diffusion 

in  the  plastic  tree  canopy . 35 


WIND  TUNNEL  STUDIES  OF  THE 
AIR  FLOW  AND  GASEOUS  PLUME  DIFFUSION 
IN  THE  LEADING  EDGE  AND  DOWNSTREAM 
REGIONS  OF  A  MODEL  FOREST 

by 

R.  N.  Meroney*  and  E.  T.  Yang** 

1.  INTRODUCTION 

Wind  movement  within  forest  stands  and  in  their  boun¬ 
dary  regions  dominates  the  exchange  processes  which  occur 
within  the  vegetative  canopy.  The  structure  of  the  timber 
stand  interacts  with  the  prevailing  winds  to  determine  fire 
spread  rates,  snow  pack,  soil  erosion,  dispersal  of  seed 
for  forest  regeneration,  blow  down,  and  rates  of  carbon 
dioxide  and  water  vapor  exchange  during  plant  metabolism. 

As  early  as  1893,  Metzger,  a  German  scientist,  in 
vestigated  the  effects  of  wind  action  on  trees.  Subse¬ 
quently,  a  variety  of  studies  have  been  made  of  the  behavior 
of  winds  well  inside  a  forest  (Bayton,  1963;  Cooper,  1965; 
Denmead,  1964;  Fons,  1940;  Huston,  1964;  Porpendiek,  1S49; 
Tiren,  1927;  Tourin  and  Shen,  1966).  Some  measurements  are 
available  for  the  variation  of  the  wind  at  the  edge  of  a 
forest  (Iizuka,  1952;  JReifsnyder,  1955).  These  measurements 
have  provided  a  rough  picture  of  a  highly  complex  and  tur¬ 
bulent  flow  field  within  the  vegetative  canopy. 

Agricultural  meteorologists,  atmospheric  scientists, 
and  many  hydrologists  are  interested  in  the  evaporation  and 
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exchange  processes  which  occur  in  vegetative  canopies.  Such 
information  permits  calculation  of  the  efficiency  of  water, 
energy,  and  C02  transport  in  plant  metabolism  and  the 
movement  of  foreign  additives  into  or  out  of  the  bulk 
of  a  canopy.  Since  1937,  experimenters  have  made  measure¬ 
ments  of  velocity,  temperature,  evaporation  rates,  and 
energy  balance  within  and  above  such  conf  gurations 
(Penman  and  Long,  I960;  Inoue,  1963;  Uchijima  and  Wright, 
1964;  Lemon,  1962).  These  measurements  have  provided  a 
rough  picture  of  a  highly  complex  and  turbulent  flow  field 
within  vegetation. 

Past  measurements  of  diffusion  from  point  or  line  source 
in  forest  configurations  seem  to  have  been  limited  to  measure 
ments  of  an  instantaneous  line  source  over  a  tropical 
rain  forest  by  Bendix  (Baynton,  1963)  ,  of  point  and  line 
source  distributions  over  a  deciduous  forest  by  Litton 
Systems  (Tourin  and  Shen,  1966) ,  of  instantaneous  point 
sources  in  a  jungle-like  deciduous  forest  by  MELPAR  (Allison, 
et  al.,  1968),  and  of  rates  of  particulate  dispersion  in  a 
forest  canopy  at  Brookhaven  (Raynor,  1967,  1969'.  These 
measurements  are  extensive  and  well  documented;  however, 
they  must  be  normalized  to  some  simplified  geometry  in  order 
to  determine  the  universal  cnaracter isti  ^  and  governing 
parameters  of  vegetative  penetration  by  a  diffusing  plume. 

Since  field  measurements  are  not  i_asy  to  obtain  because 
of  the  cost  of  providing  a  perfect  measuring  station  and 
the  difficulty  of  obtaining  cooperative  weather,  a  laboratory 
prograi  of  modeling  the  flow  i n  and  above  plant  covers  has 
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been  initiated  at  the  Fluid  Dynamics  and  Diffusion  Labora¬ 
tory  at  Colorado  State  University.  Previous  results  from 
this  program  have  been  published  by  Quarishi  and  Plate  (1965), 
Meroney  and  Cermak  (196?) ,  and  Meroney  (1968) . 

The  purpose  of  this  report  is  to  discuss  seme  measure¬ 
ments  of  diffusion  from  a  continuous  point  source  in  and 
above  a  model  forest  canopy.  The  results  of  this  stu’y  will 
consist  of : 

1)  A  description  of  the  diffusion  process  in  and  above 
the  simulated  canopy; 

2)  A  description  of  the  vertical  dispersion  of  the 
tracer  materials; 

3)  A  determination  of  the  effect  of  the  initial  fetch 
of  the  forest  canopy  on  tracer  dispersion,  and  finally, 

4)  A  determination  of  the  vertical  distribution  of  the 
eddy  diffusion  coefficients  in  and  above  the  modeled  canopy. 
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2.  MODELING  OF  A  FOREST  CANOPY 

The  wind  tunnel  has  been  used  repeatedly  by  the  forest 
meteorologist  in  his  effort  to  understand  the  complex  pattern 
of  flow  generated  by  the  tree--a  permeable,  random  shaped, 
elastic  object.  Tiren,  in  1927,  attempted  to  estimate  crown 
drag  from  conifer  branch-drag  measurements  made  in  a  wind 
tunnel  as  part  of  his  study  of  stem  form.  Wind-breaks  have 
been  studied  by  models  to  determine  soil  erosion  and  blow 
down  characteristics. 

Researchers  have  modeled  forest  behavior  using  live 
tree  boughs,  cotton  balls,  wcoden  pegs,  plastic  strips,  and 
even  wire  mesh  (Hi rata,  1953;  Irzuka,  1956;  Malina,  1941; 
Woodruff  and  Zingg,  1952).  These  studies  were  all  conducted 
to  deduce  the  qualitative  behavior  of  tree  barriers  for 
specific  problems.  The  investigators  apparently  made  no 
attempt  to  scale  dynamically  the  character  of  a  l.xve  tree 
except  to  compensate  intuitively  for  shape  and  porosity. 

To  model  completely  the  complex  geometry  and  structural 
characteristics  of  a  live  tree  is  obviously  not  practical; 
however,  measurements  made  on  coniferous  and  deciduous  trees 
in  the  wind  tunnel  and  in  the  field  suggest  that  equivalence 
of  drag  and  wake  characteristics  between  model  and  prototype 
trees  should  be  sufficient  to  study  the  general  flow  phenom¬ 
enon  (Lai,  1955;  Rayner,  1962;  Sauer  et  al.,  1951;  Walske 
and  Fraser ,  1963) . 

correlation  of  the  measurements  mentioned  above  plus 
additional  ones  made  on  live  trees  at  Colorado  State  Uni¬ 
versity  indicates  th  t  the  drag  coefficient  CD  may  vary 
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with  wind  speed  from  1.0-0. 3  (Burgy,  1961)  (Fig.  5).  These 
measurements  indicate  that  the  flow  is  inert ially  dominated 
(i.e.,  Reynolds  number  independent),  but  that  self¬ 
streamlining  of  the  tree  at  high  velocities  can  reduce  the 
effective  cross-sectional  area  for  the  more  flexible  species. 

Measurements  made  behind  small  specimens  of  Colorado 
spruce,  juniper^and  pine  trees  revealed  that  linear  wake 
growth  exists  behind  all  trees,  that  the  wake  shadows  of 
individual  branches  disappear  within  1-2  tree  crown  diameters 
downstream,  and  that  the  velocity  defect  becomes  Gaussian 
within  3-4  crown  diameters  (Fig.  6) . 

After  studying  a  variety  of  plastic,  mer.al  and  brush 
model  trees,  a  model  made  from  plastic  simulated-evergreen 
boughs  was  selected.  The  model  trees  chosen  have  an  average 
height  of  18  cm,  a  stem  height  of  5  cm,  and  a  crown  diameter 
of  7  cm.  The  model  tree  has  a  drag  coefficient  of  0.72 
over  the  velocity  range  studied  and  a  lateral  wake  growth 
similar  to  that  measured  for  live  trees  (Figs,  5  and  6). 

Results  of  extensive  single  tree  drag  measurements  made 
within  regular  geometric  arrays  of  the  same  model  tree  (an 
orchard  arrangement)  are  reported  by  Hsi  and  Nath  (1968) . 

The  drag  p,. ''files  measured  show  a  similar  behavior  to  the 
bending  moment  measurements  made  by  Walske  and  Fraser  (1963); 
that  is,  there  is  a  sharp  decrease  in  drag  on  the  trees  with 
distance  down-wind  followed  by  a  slight  rise  to  an  asymp¬ 
totic  constant  value. 

Shear  plate  measurements  made  within  the  random  canopy 
array  under  discussion  herein  display  the  same  characteristics 
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as  the  regular  arrangements.  Figure  (7)  plots  local  shear 
force  vs  distance  downwind  from  the  canopy  inception.  The 
minimum  observed  within  the  first  2  m  is  evidently  the  re¬ 
sult  of  a  relatively  stagnant  region  inside  the  canopy  which 
also  explains  the  behavior  of  the  diffusion  plume  discussed 
subsequently.  This  same  phenomenon  was  found  for  flow  over 
a  mouel  peg  canopy  (Merouey  and  Cermak,  1967) . 
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3.  EXPERIMENTAL  EQUIPMENT  AND  PROCEDURES 

3.1  Wind  Tunnel  and  Canopy  Arrangement;  The  experimental 
data  were  obtained  in  the  low  speed  Army  Meteorological  Wind 
Tunnel  in  the  Fluid  Dynamics  and  Diffusion  Laboratory  at 
Colorado  State  University  (Plate  and  Cermak,  1963).  This 
tunnel  was  specifically  designed  to  study  fluid  phenomena  of 
the  atmosphere.  The  tunnel  has  a  2  m  square  by  26  m  long 
test  section  with  an  adjustable  ceiling  to  provide  a  zero 
pressure  gradient  over  the  forest  canopy.  The  model  trees 
were  inserted  into  holes  in  aluminum  plate  sections  which 
extended  the  width  of  the  tunnel  and  11  m  downstream  from 
the  tunnel  midsection.  The  elements  were  randomly  position¬ 
ed  with  approximately  one  tree  per  36  cm2.  From  above, 
this  arrangement  gave  the  same  visual  appearance  as  a 
moderately  dense  coniferous  forest.  This  density  would  be 
equivalent  to  a  stand  density  index  as  calculated  by  Reinke 
(1933)  of  250  for  a  forest  with  an  average  tree  height  of 

40  ft  and  a  diameter  at  breast  height  of  10  inches  (Fig.  1), 
(Fig.  2) .  A  volumetric  density  number  has  been  calculated 
to  describe  the  canopy  density  by  Sadeh,  et  al.,  (1969). 

When  one  describes  the  volume  occupied  by  a  single  tree  as 
a  combination  of  a  crowncone  and  trunk  cylinder,  the  ratio 
of  tree  occupied  volume  to  volume  beneath  the  mean  canopy 
height  is  26%. 

3.2  Velocity  and  Turbulence  Measurements :  A  single  wire 
constant  temperature  anemometer  was  used  to  measure  velocity, 
turbulent  intensity,  and  shear.  In  addition,  pitot-static 
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tube  measurements  were  made  at  each  section.  The  sensing 
elements  of  the  anemometer  circuit  were  platinum  wire  0.2 
mil  in  diameter  and  approximately  0.25  cm  long.  The  bridge 
circuit  utilized  was  a  CSU  Solid  State  Anemometer.  The 
pitot  tube  output  went  to  a  Transonic  Model  A,  Type  120 
electronic  pressure  meter.  Turbulence  signals  were  inter¬ 
preted  by  means  of  a  Bruel  and  Kjaer  RMS  meter,  Model  2416. 

3 . 3  Concentration  Measurement-Helium  Tracer  Gas:  The 
character  of  the  flow  field  was  studied  by  mapping  the  dif¬ 
fusion  plume  of  a  continuous  point  source.  Helium  gas  was 
used  as  one  tracer  for  the  diffusion  experiment.  The  gas 
was  released  continuously  at  a  constant  rate  of  630  cc/min 
from  a  2  mm  nozzle  located  in  or  above  the  canopy.  The 
sampling  probe,  manufactured  from  small  diameter  hypodermic 
tubing,  was  mounted  on  a  traversing  carriage,  the  horizontal 
and  vertical  positions  of  which  were  controlled  remotely 
from  outside  the  tunnel.  Helium  concentration  was  measured 
at  ground  level  along  a  line  normal  to  the  axis  of  the  plume 
and  vertically  at  the  plume  centerline. 

Samples  were  dtown  into  the  probe  at  a  constant  rate 
and  passed  over  a  standard  leak  into  a  mass  spectrometer 
(Model  of  the  Vacuum  Electronic  Corporation)  .  Output 

of  the  mass  spectrometer  was  an  electrical  voltage  pro¬ 
portional  to  concentration.  The  mass  spectrometer  was  cal¬ 
ibrated  periodically  by  a  set  of  pre-mixed  gases  of  research 
grade.  Figure  3  shows  the  experimental  arrangement. 
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Since  a  closed-circuit  wind  tunnel  was  used,  the  ambient 
concentration  level  of  helium  built  up  in  the  wind  tunnel  with 
time.  Eventually,  most  of  the  gas  did  leak  out;  therefore  the 
amount  of  helium  in  tne  ambient  flow  was  never  higher  than  60 
parts  per  million.  Nevertheless,  an  ambient  concentration 
measurement  was  taken  after  each  profile.  The  relative  con¬ 
centration  was  obtained  by  subtracting  the  corresponding  am¬ 
bient  concentration  from  the  absolute  concentration.  All  data 
presented  in  the  figures  or  tables  are  relative  concentrations. 

Due  to  the  slow  response  of  the  mass  spectrometer,  a  per¬ 
iod  of  one  to  two  minutes  was  allocated  for  the  stabilization 
of  each  reading  before  it  was  recorded.  Usually,  the  con¬ 
centration  signal  itself  was  averaged  over  at  least  60  seconds. 
This  method  gave  results  that  compared  favorable  with  the 
average  of  signals  taken  over  a  period  as  long  as  250  seconds 
by  graphical  means. 

3 . 4  Concentration  Measurement  -  Kr-85  Tracer  Gas:  To 
investigate  the  buoyancy  character  of  the  helium  -racer 
additional  measurements  were  obtained  utilizing  a  mixture 
of  Kr-85  and  air  as  a  tracer.  It  is  a  radioactive  noble  gas 
which  does  not  chemically  combine  with  any  other  molecules 
in  the  system  studied.  Krypton-85  has  a  half  life  of  10.6 
years  so  there  is  no  appreciable  decay  during  a  diffusion 
experiment.  The  radioactive  gas  was  diluted  about  a  million 
times  before  use  and,  as  such,  has  physical  characteristics 
equivalent  to  those  of  air .  Its  detection  procedure  is 
fairly  simple  and  direct.  Handling  and  safety  procedures 
for  wind  tunnel  experiments  with  Kr-85  tracer  gas  have  been 
discussed  in  detail  by  Chaudhry  and  Meroncy  (1969). 
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The  flow  rate  of  Kr-85  mixture  was  controlled  by  a 
pressure  regulator  at  the  bottle  outlet  and  monitored  by  a 
Fisher  and  Porter  flowmeter.  Source  concentration  was  6.4 
y-curie/cc  of  Kr-85,  a  beta  emitter. 

A  sampling  rake  of  eight  probes  was  manufactured  from 
2  mm  diameter  hypodermic  tubing  and  was  mounted  on  a 
traversing  carriage  whose  horizontal  and  vertical  position 
was  controlled  remotely  from  outside  the  tunnel.  Concentra¬ 
tions  were  measured  at  ground  levels  at  various  scaled  dis¬ 
tances  from  200  to  400  feet  downwind  and  at  vertical  eleva¬ 
tions  centered  on  plume  maximum  concentrations.  Samples  were 
aspirated  at  a  constant  rate  of  500  cc/min  into  eight  TGC-308 
Tracerlab  Geiger-Mueller  side  wall  cylindrical  counters. 

Samples  were  flushed  through  the  counting  tubes  for  at  least 
two  minutes,  Valve  A  in  pigure  (5B)  was  closed,  and  each  sample 
was  subsequently  counted  for  one  minute  on  Nuclear  Chicago 
Ultra-scaler  Model  192A.  All  samples  counted  were  adjusted 
for  background  radiation  (See  Fig.  4a  and  4b) . 
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4 .  EXPERIMENTAL  RESULTS 

All  measurements  were  taken  at  a  free  stream  velocity 
of  6  ra  sec  ^ .  The  ceiling  of  the  test  section  was  adjusted 
for  zero  pressure  gradient  and  the  upstream  velocity  profile 
was  measured  and  found  to  be  logarithmic.  The  temperature 
condition  was  constant  and  hence  neutral  stability  existed. 

4 . 1  Typical  Velocity  and  Turbulent  Intensity  Profile  Results; 
A  sequence  of  vertical  profiles  of  mean  velocity  measurements 
were  made  along  the  tunnel  centerline  both  in  and  above  the 
forest  canopy.  The  transformation  of  the  wind  profiles  in 
the  vertical  direction  are  shown  in  Figure  (5) .  Jetting  of 
the  wind  flow  beneath  the  canopy  is  observed  for  at  least  the 
first  3  m  (or  15  canopy  heights);  subsequently,  the  wind 
profile  reaches  an  equilibrium  state  at  about  4  m  (or  20 
canopy  heights) .  Finally,  accelerations  of  the  wind  are 
observed  during  the  last  2  m  of  the  canopy  as  the  wind  ad¬ 
justs  to  the  smooth  surface  downwind.  The  extent  of  the 
entrance  region  agrees  with  previous  measurements  by  Meroney 
and  Cermak,  and  Plate  and  Quarishi  (1965),  but  is  greater 
than  that  tentatively  suggested  by  Reifsnyder  (1955) .  The 
shape  of  the  equilibrium  velocity  profile  agrees  qualitative¬ 
ly  with  prototype  measurements  for  moderately  dense  conifer 
forests  (Cooper,  1965;  Denmead,  1964;  Fons,  1940;  Poppendiek, 
1949;  Reifsnyder,  1955;  Tiren,  1927;  Tourin  and  Shen ,  1966). 

In  the  winter  the  Minnesota  deciduous  forest  of  Tourin 
and  Shen  (1966) ,  compares  favorably  quantitatively  with  a 
fairly  dense  peg  arrangement  (Fig.  10),  whereas,  the  plastic 
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tree  canopy  simulates  summer  measurements  made  by  Allen 
(1968),  Shinn  (1969),  and  Tourin  and  Shen  (1966),  (F.ig.  11). 

Velocity  data  from  the  plastic  tree  canopy  has  also 
been  compared  with  prototype  measurements  by  means  of  a 
dimensionless  velocity  defect  argument.  Shinn  (1969)  cal¬ 
culated  the  defect  between  the  pre-canopy  velocity  profiles 
and  that  measured  within  the  forest.  The  result  for  a  fetch 
length  of  x/h  =5  is  displayed  in  Figure  (12). 

The  profiles  above  the  canopy  are  logarithmic  and  can 

be  plotted  to  follow  the  displacement  law  u/u*  =  k_1ln[(y-d) 

/zQ)  as  shown  by  Plate  and  Quarishi  (1965)  .  However,  it 

should  be  noted  that  the  popular  regression  technique  first 

suggested  by  Lettau  to  solve  for  u*,  d,  and  2  could  not 

o 

be  utilized  unless  modified  (Robinson,  1961).  This  program 
(a  version  of  which  is  known  as  the  "Three  Bears"  program)  un¬ 
fortunately  assumes  u*,  d,  and  zq  are  independent;  as  a 
result,  some  investigators  have  obtained  the  physically 
suspect  result  that  d  is  negative  (Rung,  1961).  In  our 
computations,  d  was  assumed  equal  to  the  canopy  height; 
thus  zq  -  22  cm,  and  u*  =  14  m/sec.  In  addition,  measure¬ 
ments  over  the  peg  canopy  suggested  that  the  velocity  pro¬ 
files  may  be  dominated  by  the  canopy  top  wake  until  z  -  2.5 
to  3  h;  hence,  it  would  appear  that  forest  micro-meteorolo¬ 
gists  should  not  attempt  a  log-law  analysis  unless  they 
utilize  fairly  tall  towers.  Moreover,  recent  analysis  of 
data  for  above  canopy  flows  suggests  that  tho  friction 
velocity  and  roughness  length  are  not  local  quantities  but 


vary  with  height;  perhaps  because  the  assumption  of  a  constant 
shear  stress  region  is  invalid,  (Sadeh,  et  al.,  1969). 

Hot  wire  anemometers  were  used  to  measure  turbulence  char¬ 
acteristics  in  and  over  the  model  canopy  (Fig.  9).  Values  of 
longitudinal  intensity  up  to  0.35  were  measured  in  and  above  the 
model  forest  canopy.  They  correspond  to  field  measurements  by 
Tourin  and  Shen  (1966)  who  report  average  values  of  longitudinal 
turbulence  of  0.33  at  the  40  foot  level.  Subsequent  measurements 
by  (Sadeh,  et  al.,  1969)  also  measured  high  turbulence  intensity 
levels;  however,  changes  in  measurement  techniques  resulted  in 
values  as  high  as  0.77  in  the  established  flow  regime.  Tourin 
and  Shen  also  noted  the  decrease  of  turbulence  as  one  moves 
downward  into  the  forest  cover. 

4.2  Diffusion  Plume  Results:  Plumes  were  released  at  the  model 
forest  entrance  from  locations  near  the  ground,  at  half  canopy 
height,  and  at  the  top  of  the  canopy.  Releases  were  also  made 
in  the  equilibrium  wind  profile  region  downstream.  Tables  1 
through  7  summarize  data  measured. 

Figures  (13)  and  (14)  display  the  typical  plume  exhalation 
by  the  forest  near  the  entrance  and  the  subsequent  re-inhalation 
further  downstream.  A  similar  behavior  has  been  noticed  for  re¬ 
leases  of  gas  over  a  model  crop  canopy  simulated  with  dowel  pegs 
(Meroney  and  Cermak,  1967,  Yano,  1967).  This  phenomena  is  a  re¬ 
sult  of  vertical  motions  near  the  front  of  the  forest  canopy  pre¬ 
viously  reported  by  Iizuka  (1S52).  The  subsequent  rapid  penetra¬ 
tion  further  downstream  may  be  due  to  the  intense  shear  and  mix¬ 
ing  near  the  canopy  top  over  the  initial  fetch  region.  The 
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ramification  of  this  effect  upon  fire  spread  and  parasite  con¬ 
trol  by  spray  is  obvious. 

Plume  releases  within  the  forest  near  the  ground  were  char¬ 
acterized  by  wide  meandering  and  large  lateral  dispersal.  Such 
erratic  behavior  including  plume  bifurcation  occurs  frequently 
during  forest  diffusion  experiments  (Allison,  1968;  Shinn,  1969; 
Geiger,  1950). 

Figures  (15)  through  (18)  present  vertical-isoconcentration 
sections  through  continuous  point  source  plumes  released  at  var¬ 
ious  heights  above  the  ground  (i.e.,  0,  l/2h,  h,  and  l-l/2h) 
where  the  flow  field  appears  fully  established  (i.e.,  x/h  =  33). 
For  the  elevated  releases  the  sequence  of  stages  of  the  concen¬ 
tration  gradient  observed  upon  penetration  of  the  plume  down¬ 
stream  are  similar  to  those  observed  by  Flemming  (1967)  during 
elevated  line  source  releases  over  a  deciduous  forest.  Initially, 
there  is  a  qrad^nt  downward  followed  by  a  gradient  in  concen¬ 
tration  upward  even  farther  downstream. 

It  is  interesting  to  note  how  the  diffusing  cloud  tilts 
forward  near  the  tree  top  due  to  wind  shear,  and  how  a  rapid 
forward  movement  has  resulted  from  the  relatively  high  wind 
speed  at  the  tree  tops.  The  very  rapid  vertical  growth  of  the 
plume  for  ground  source  releases  is  another  feature  also  dupli¬ 
cated  by  ground  based  boniblet  measurements  (Tourin  and  Shen, 

1969).  Tne  MELPAR  study  did  not  incorporate  ary  significant 
number  of  vertical  measurements;  however,  observation  of  putf 
behavior  led  to  the  conclusion  vertical  mixing  to  the  canopy 
top  was  complete  within  very  short  donwind  distances  (MELPAR , 


1968)  . 
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It  has  been  generally  observed  for  continuous  plume  re¬ 
leases  that  the  maximum  concentration  at  ground  level  decreases 
at  a  rate  proportional  to  a  power  function  of  the  longitudinal 
downstream  distance,  x~ra.  For  a  plume  dispersing  in  or  above 
a  vegetative  canopy,  the  rate  of  dispersal  also  appears  to  be 
a  function  of  the  distance  from  the  release  position,  (x-xs)  m, 
(see  Fig.  19).  The  rate  of  dispersion,  however,  is  much  larger 
than  for  plumes  dispersing  over  a  smooth  surface  (Malholtra 

and  Cermak,  1964),  (i.e.,  mcanopy  =  -4.8,  raper  *  “2.5, 

plastic  canopy 

m  —  _  i  c  \ 

smooth  • D} • 

surface 

Examination  of  bomblet  releases  in  a  deciduous  forest  by 
Tourin  and  Shen  (1966)  produced  values  of  m  =  -7.0  for  a 
typical  near-neutral  summer  release  and  m  =  -3.0  for  a 
winter  release.  The  average  decay  rate  for  all  F.P.  releases 
in  a  summer  jungle  canopy  was  found  to  be  -3.1  by  MELPAR, 

Inc.  (1968). 

Brown,  et  al .  (1969)  have  proposed  that  the  ventilation 
rate  of  most  vegetative  canopies  may  be  correlated  to  an  en¬ 
vironmental  index  defined  as  El  =  u  /  u,  where 

ac.  be. 

u  ,  =  velocity  at  two  canopy  heights. 

=  velocity  at  one-half  canopy  heights. 

If  the  coefficient  -m  is  plotted  versus  such  an  environ¬ 
mental  index  one  notes  an  increase  in  dispersion  rate  a:;  the 
index  increases  followed  by  a  decrease  tc  zero  for  very  dense 
vegetative  configurations.  Th ' s  behavior  appears  to  cor¬ 
relate  with  the  increase,  m  turbulent  intensity  initially 
until  tne  blockage  becomes  so  great  as  to  inhibit  the  rate 
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of  disperson  of  the  gases,  after  which  -m  decreases,  see 
Figure  (26). 

When  the  flow  above  and  below  the  canopy  ceiling  are 
treated  as  separate  flow  regimes,  similarity  conditions 
appear  to  exist  when  the  appropriate  characteristic  length 
parameters  are  chosen.  If  the  character  of  the  concentration 
profile  is  examined  above  the  canopy  top,  one  finds  that 
similarity  may  be  obtained  over  long  fetch  distances  by  dis¬ 
playing  C/C^  vs  (z  -  h)/(A  -  h) ;  where  h  =  canopy  height, 
and  \  =  characteristic  width  of  plume  when  C  =  %  C^, 

(Fig.  20) .  Data  is  compared  to  an  analytic  expression  which 
also  summarize-  the  character  of  plume  releases  over  smooth 
surfaces. 

Comparison  of  isoconcentration  profiles  for  the  Helium 
tracer  gas  and  Kr-85  tracer  gas  suggests  that  the 
initial  buoyancy  of  the  undiluted  Helium  source  had  little 
effect  on  the  dispersion  in  and  above  the  canopy.  Figures 
(16a)  and  (16b)  display  the  measurements  for  the  Helium  and 
Krypton  tracers  respectively.  In  addition,  slight  variations 
observed  in  the  ground  level  concentration  variation  with 
downward  distance  are  not  of  the  order  or  directicr.  to  be 
attributed  to  buoyancy  effects. 

4.3  Cddy  Diffusion  Coefficient:  The  concept  of  a  macro¬ 
scopic  equation  of  turbulent  dispersion  of  some  property  C 
results  generally  in  the  equation 
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where  Kv  is  the  coefficient  of  turbulent  diffusion.  The 
xi 

coefficient  Kx,  incorporates  within  itself  the  . omplexities 
of  the  actual  transport  process.  Hence,  most  analytical 
studies  of  fluid  mechanics  require  some  theo-etical  or  em¬ 
pirical  expression  for  the  variation  of  Kx^  with  other 
parameters.  Several  scientists  have  studied  the  nature  of 
Kx^  for  plant  communities,  but  further  data  are  still  needed 
(Penman  and  Long,  1969;  Inoue,  1963;  Yano,  1966;  Saito,  1964). 

The  oddy  diffusion  coefficient  for  transport  of  the  in¬ 
jected  gas  in  the  model  canopy  has  been  determined  utilizing 
concentration  and  velocity  profiles  and  a  finite  difference 
interpretation  of  Equation  (1).  In  order  to  simplify  the 
discretization  analysis  the  concentration  data  were  converted 
to  line  source  data  by  the  assumption  of  normal  distributions 
and  lateral  integration.  Two  computational  methods  were 
utilized  to  calculate  Kz(z).  In  one,  Equation  (1)  was 
solved  directly  in  finite  difference  form  for  K2(z)  such 
that 
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(1)  was  integrated  once  in 

rivative  term  such  that 


Kz(z)  = 


are  replaced  by  their  finite 
In  the  second  method,  Equation 
z  to  eliminate  the  second  de- 
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These  methods  gave  essentially  identical  results  in  and  above 

the  forest  canopy.  Calculations  were  performed  on  a  CDC  6400 

computer  at  Colorado  State  University  using  input  data  taken 

from  lines  faired  through  the  ground  source  concentration 

measurements,  at  x  =  6  Meters  and  from  vertical  velocities 

s 

calculated  from  the  slope  of  streamlines. 

The  resulting  profiles  in  K{z)  are  displayed  in  Figure 
(21).  Three  distinct  regions  of  variation  of  K  are  notice¬ 
able.  Immediately  adjacent  to  the  wall  is  a  zone  where  K 
increases  exponentially.  In  the  area  from  4  to  12  cm,  K 
remains  essentially  constant;  and  K  becomes  proportional 
to  (z-d)  where  d  is  a  displacement  height.  Similar  be¬ 
havior  has  been  observed  for  prototype  canopies.  Finally, 
these  K  profiles  may  also  be  described  as  qualitatively 
similar  to  the  peg  data. 

A  number  of  authors  have  suggested  that  K  should  re¬ 
main  constant  in  vegetative  cover;  others  have  suggested 
that  K  should  vary  linearly  (Inoue,  1963;  Uchijima  and 
Wright,  1964).  It  is  interesting  to  note  that  for  the  case 
of  the  model  peg  canopy,  both  conditions  of  K  exist,  al¬ 
though  in  different  regions.  Figure  (22)  comp?«-*'S  the  dis¬ 
tribution  of  K  within  the  canopy  with  typical  results  of 
the  distribution  of  K  for  a  pine  forest  as  measured  by 
Denmean,  (1964). 

The  experimental  data  mesh  from  whi~b  the  estimates  of 
K^fz)  were  obtained  was  fairly  coarse;  hence,  to  verify 
the  results  it  was  decided  to  recompute  the  concentration 
distributions  numerically  for  the  elevated  release  conditions 
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for  a  continuous  point  source  situation.  Equation  (1)  was 

discretized  and  solved  by  means  of  an  alternating-direction- 

implicit  technique  described  by  Peacestan  and  Rachford  (1955) . 

Initially  it  was  assumed  K  =  K  (z) . 

y  z 

Figure  (23)  compares  the  ground  concentrations  as 
measured  and  as  calculated  when  initial  plume  concentrations 
at  x  ■  25  cm  were  substituted  into  the  calculation  proce¬ 
dure.  If  a  value  of  the  ratio  K y/Kz  =  2.0  or  4.0  is 
assumed,  one  obtains  a  somewhat  better  comparison  as  shown 
on  the  same  figure.  The  value  of  is  normally  expected 

to  exceed  K  especially  in  the  near  ground  region.  Faster 

Z 

lateral  dispersion  at  ground  level  has  also  been  observed 
for  model  peg  canopies  (Meroney  and  Cermak,  1967). 

Figure  (24)  displays  the  result  of  the  assumption 
K y/Rz  1  1  upon  the  cross-section  isoconcentrations  lines  as 
seen  for  an  elevated  and  ground  release  in  the  plastic  tree 
canopy . 

4.4  Forest  Penetration  Model:  Despite  the  existence  of  com¬ 
plex  sets  of  diffusion  data  in  various  vegetative  canopy  con¬ 
figurations,  only  elementary  solutions  for  understanding 
physical  dispersion  of  gases  in  -orests  has  been  put  forward. 
Most  experimentalists  have  tried  to  fit  their  results  to  re¬ 
gression  equations  (Baynton,  1963;  Tourin  and  Shen,  1969; 
Allison,  1968);  for  example  Baynton  (1963)  suggested 

(DO*^' 'ground  ■  IA  *  ~C--5u--EVf 1  tDos^el above  canopy, 

where  U  is  velocity  above  the  canopy,  ,‘.T  is  temperature 
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difference  above  and  below  canopy,  and  o.  is  standard 
deviation  of  wind  direction  above  forest.  As  Baynton  notes 
such  a  formula  applies  specifically  to  the  forest  in  which  th«5 
data  were  collected  since  the  height  of  the  forest  and  forest 
density  are  not  parameters.  Baynton  could  detect  no  below 
canopy  mean  and  speed  in  his  dense  jungle  canopy;  hence  his 
regression  formula  only  allows  for  vertical  diffusion  in 
and  out  of  the  forest  with  no  longitudinal  convection. 

Tour in  and  Shen,  on  the  other  hand,  worked  in  a  somewhat  less 
dense  canopy  and  suggested  that  the  relation 

(Dosage)  ground  =  0  51  x“0.993  -  -0.75  -  -0 . 98 (]_F) 0. 25 

where  °c  s  standard  derivation  of  vertical  angle  at  the 
40  meter  level, 
u  =  mean  and  speed,  and 

F  =  tree  canopy  density  based  on  light  intensity 
measurements  yielded  the  best  fit  to  all  available  line 
source  data.  The  longitudinal  decay  parameter  from  the 
Litton  Systems  study  of  -0.993  compares  with  a  value  of  -0. 
for  this  work.  in  addition  to  modifications  of  simple 
Gaussian  plume  models  (Tourin  and  Shen,  1969)  (Allison,  1968), 
one  may  also  appeal  to  a  simple-minded  one-dimensional  model 
for  canopy  penetration,  first  suggested  by  Calder,  (1961). 

The  below  canopy  concentrations  resulting  from  an 
elevated  continuous  release  line  source  can  be  estimated  by. 
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where  s  *  penetration  coefficient  and  u  =  below  canopy 
wind  speed.  The  above  canopy  measurements  have  been  fitted 
to  the  formula  suggested  by  Bosanquet  and  Pearson  (1936), 

Cabove  (x)  *  £exp(-  B/x) , 
canopy 

and  the  predicted  below  canopy  concentrations  compared  with 
experimental  data  in  Figure  (25) .  Obviously  the  Bosanquet 
formula  is  somewhat  inadequate,  however,  it  is  apparent  fair 
comparison  is  obtained  for  a  model  penetration  coefficient 
of  0.75  sec  This  is  comparable  to  a  prototype  exchange 
rate  of  '0.45  minutes  1  since  the  time  scale  for  the  model 
may  be  interpreted  as  100  times  less  than  in  the  field. 

Calder  also  suggested  a  manner  in  which  to  check  the 
validity  of  the  mathematical  model  and  estimate  the  parameter 
H  *  s/u.  He  noted  that  the  model  requires  that 
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for  different  selected  values  of  the  transform  parameter  p. 
This  equation  was  checked  numerically  for  a  range  of  p 
from  2  to  10,  and  die  calculated  parameter  H  varied  from 
1.92  to  1.14;  whereas,  the  best  first  value  from  the  figure 
appears  to  be  1.50. 

Although  the  model  for  an  instantaneous  point  source 
suggested  by  the  MELPAR  (i9b3)  study  incorporated  vertical 
and  lateral  dispersion  degrees  of  freedom  tneir  predictions 
were  limited  to  below  canopy  release  conditions.  Ir.  addition. 
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they  incorporated  an  infinite  mass  sink  at  the  canopy  top, 
which  was  admitted  to  be  over  restrictive.  Information  con¬ 
cerning  the  vertical  concentration  profiles  obtained  in  this 
study  might  be  used  to  improve  the  MELPAR  model,  since  no 
vertical  measurements  were  available  in  the  Jungle  Canopy 
study. 

Tourin  and  Shen  also  compared  their  measurements  for 
elevated  line  source  releases  above  a  Wisconsin  forest  with 
Calder's  model  and  another  model  developed  from  Lattau's 
hypothesis  of  vorticity  transfer.  These  models  generally 
did  not  agree  with  the  observed  data,  as  well  as  the  regression 
equation;  however,  one  can  not  tell  whether  this  is  a  failure 
of  the  below  canopy  models  utilized  or  the  inadequacy  of  the 
Bosanquet-Pearson  expression  used  to  predict  above  canopy 
dosages. 
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5.  CONCLUSIONS 

It  is  apparent  that  the  general  character  of  flow  in 
and  above  vegetative  canopies  may  be  satisfactorily  simulated 
in  the  meteorological  wind  tunnel.  In  addition,  these  new 
data  suggest  that  even  the  micro-structure  transport  phenomena 
behave  in  a  manner  similar  to  that  of  the  prototype.  There¬ 
fore,  it  is  possible  to  conclude  that: 

1)  The  basic  trends  of  the  dynamic  and  kinematic  be¬ 
havior  of  a  complex  vegetative  cover  may  be  simulated  by  a 
simple  porous  geometry  in  a  wind  tunnel. 

2)  The  initial  fetch  of  the  peg  canopy  affects  tracer 
dispersion  of  a  continuous  point  source  in  a  unique  manner: 
Vertical  convective  motions  exhale  the  gases  released  at 
the  beginning  of  the  canopy,  and  subsequently,  the  canopy 
appears  to  re-inhale  the  products  farther  downstream. 

3)  The  concentration  profile  above  the  canopy  displays 
the  features  of  a  plume  released  over  a  flat  plate  but  dis¬ 
placed  by  a  height  h. 

4)  The  eddy  diffusion  coefficient  varies  linearly  as 
(z-d)  above  a  vegetative  cover  and  has  a  growth  rate  nearly 
propc. Lional  to  ku*. 

5)  The  eddv  diffusion  coefficient,  K  ,  within  the 

'  Z 

artificial  vegetative  cover,  appears  to  develop  into  three 
regions:  Initially  K grows  exponentially,  next  it  re¬ 
mains  constant,  and  finally,  grows  at  a  linear  rate. 

6)  The  experimental  law  for  attenuation  of  boundary 

-4 . 8 
x 


concentration  was  obtained  as 


for  gas  source  t  -'leases 
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far  from  the  canopy  inception.  (Rates  of  dispersion  are  some¬ 
what  larger  near  the  edge  of  the  vegetative  cover . ) 

7)  The  lateral  eddy  diffusion  coefficient,  ,  appears 
to  be  2  times  larger  than  the  vertical  transport  rate  as 

on  approximation.  However,  it  is  expected  that  +  0  at 
ground  level. 

8)  Considering  the  similarity  of  plume  behavior  when 
considered  separately  above  and  below  the  top  of  the  canopy, 
it  would  appear  that  models  directed  to  treat  the  physics 

of  these  two  layers  separately  are  justified. 
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Concentration  Profiles  of  Diffusion 
in  the  Plaatic  Tree  Canopy 
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figure  2.  Model  pi^ntic  forest 
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Figure  3.  Helium  detec  :ion  system. 
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Figure  5.  Drag  coefficient  of  live  and  model  trees. 


Figure  6 .  Wake  charac 
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Figure  8.  Velocity  profiles. 
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Figure  10.  Comparisons  with  winter  forests. 
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Figure  11.  Comparisons  with  summer  forests 
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Figure  12.  Velocity  defect  comparison. 


Figure  13.  Diffusion  -  Isoconcentrati^n  profiles, 
zs  =  0.0  cm 
xs  =  0.0  m 


Figure  1-.  Diffusion  -  Isoconcentration  profiles, 
zs  -  10.0  cm 
x  s  -  0.0  m 


Figure  16a.  Diffusion  -  Isoconcentration  profiles 
zs  =  10.0  cm 
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Figure  17.  Diffusion  -  Isoconcentration  profiles, 
zs  =  18.0  cm 
xs  =  6.0  m 
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Figure  21.  Eddy  di 'fusion  coefficient  -  mass. 
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figure  -4 .  Cron j  -  section- isoconcentration  profil«». 
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